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Extending the method successful for one-loop integrals, the computation 
of two-loop diagrams with general internal masses is discussed. For the 
two-loop vertex of non-planar type, as an example, we show a calculation 
related to Z° — > ti vertex. 

The calculation of loop integrals is essential to obtain the precise 
theoretical prediction for high-energy reactions. Loop integrals for one- 
loop diagrams can be expressed by logarithms and dilogarithmscl. For 
a class of higher order diagrams, compact analytic expressions are ob- 
tained. For instance, some tworloop two point functions with single 
mass(m) are given by functionsH of x = s/m 2 . As two-loop diagrams 
in the electro-weak theory include complicated mass combinations, an 
analytic formula seems not available. In general, for the case with 
more than two independent dimensionless variables, it seems to be im- 
possible to obtain a compact analytic formula by polylogarithms and 
so forth. In scattering processes, we have two or more invariants and 
in the electroweak theory we encounter diagrams with many different 
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masses. Thus the theoretical prediction in the electroweak theory re- 
quires a method to compute the loop integrals for arbitrary mass scales. 
To meet this request, we have already developed numerical methods for 
two-loop integrals. There exist a few other approaches for this direc- 
tion by momentum space integral,^ momentum expansion with Pade 
approximation^ asymptotic expansion^, and so forth. 

Automatic generation of two-loop diagrams is possible in the stan- 
dard modelQ. From the generated diagrams, symbolic expression for the 
numerator and the loop integrand can be generated. The latter, the 
integrand in two-loop as a function of Feynman parameter, is treated 
numerically 

Our strategy is as followsQ: 

• Any loop integral can be given by the definite integral in Feynman 
parameter space. 

• The integrand can be singular where the denominator becomes 
zero. If there is no singularity, the value of integral can be ob- 
tained without any trouble. The numerical accuracy is only lim- 
ited by available computer time. 

• If there exists singularity, the numerical integration is done after 
either of the following: 

1. (Symmetric method) The integrand is replaced by a sym- 
metric sum of the function so as to eliminate the singularity. 

2. (Hybrid method) The integral for a few of variables is done 
analytically until the singularity is replaced by an integrable 
function. 

In this report, we present results for those integrals appearing in the 
non-planar vertex functions. 

We consider the following scalar integral for vertex function 



where l/D^s are propagator functions. Introduction of Feynman pa- 
rameters transforms it in the following form. 
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where U and V can be determined from the topology of the diagram 
H If we confine ourselves to a class of diagrams where there is no 



Numerical Approach to Two-Loop Integrals with Masses 3 



ultraviolet divergence in a subgraph, we can take n = 4. The function 
Uis 

tf=En*i ( 3 ) 

where the summation is taken over co-trees(T) in the diagram. The 
function V is given by 



1 



V = y Ex j m>--'£W g p t g 



(4) 



where the second sum is taken over cutsets(S) in the diagram. 

External momenta, pi,P2,P3 are defined to flow inward into the ver- 
tex and they satisfy p\ + p 2 + P3 = 0. We use the notation = pj, 

Xi„, = 1 — X{.... In the following, we consider 

the integral 



where 



and 



r 1 1 
J= dxi ■ ■ ■ dxaSil — x)-t— — 
Jo v ^ '{V + it) 



XjTTlj 



w sPs = h s i + hs2 + /3S3 • 
For the diagrams in Fig.|T], we obtain the following results. 



(5) 
(6) 
(7) 





P2 Pi 

Figure 1: Two-loop vertex diagrams, (a)planar type, (b) non-planar 
type. 



(a) Planar type 



U — £12X3456 + £3^456 



(8) 
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fl = X 123 X4X 6 + X1X3X5 

h = + x 2 x 3 x A (9) 

h = XlX 2 X U 56 + XU3X4X5 + XiX 3 X 5 + X2X3X4 

(b) Non-planar type 

U = X12X3456 + x u x 56 (10) 



fi = X 1256 X 3 X 4 + XiX 3 X 6 + X 2 X 4 X 5 

h = X1234X5X6 + x 2 x 3 x 6 + X1X4X5 (11) 

fi = X 34 56XlX 2 + X 2 X 4 Xq + X1X3X5 

For the first case, we have already reported the nnmerical results!. 
Here, we report the results for the non-planar vertex. 
First, variables are transformed as follows. 

xi = z 3 (l + y 3 )/2, x 2 = z 3 {\ -y 3 )/2 

x 3 = z 1 {l-y 1 )/2, X4 = z 1 {l+y 1 )/2 (12) 

x 5 = z 2 (l - y 2 )/2, x 6 = z 2 (l +y 2 )/2 

After the transform, three-fold symmetry in the non-planar vertex can 
be seen clearly. The integral becomes 

1 r 1 r 1 r 1 r 1 1 

J=- dz 1 dz 2 dz 3 S(l-Y,z)z / dy 1 dy 2 dy 3 , (13) 
8 Jo J-i J-i J -1 {D + ie) z 

where z = z\z 2 z 3 , and the function in the denominator becomes 

V = t yAy + b-y + c (14) 

where 




(15) 

V 2/3 J 

I / -z\z\8\ z(-si - s 2 + s 3 )/2 z(-si + s 2 - s 3 )/2 

A = -\ z(-si - s 2 + s 3 )/2 -z\z 2 s 2 z{s x - s 2 - s 3 )/2 

z{s x + s 2 - s 3 )/2 z(s x - s 2 - s 3 )/2 -23*3*3 

(16) 

, / z 1 (-ml + ml) 
b=-\ z 2 (-mj + ml) |, (17) 
\ z 3 {m\-m\) 



4 



c=-U 
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Z\Si + z 2 s 2 + z 3 s 3 - 2(m\ + ml)zi - 2(m| + m 2 6 )z 2 - 2(m\ + ml)z 3 

(18) 

U = z x z 2 + z 2 z 3 + z 3 z x . (19) 



So the problem is transformed into how to carry out the double 
integral of a function similar to that which appears in a box diagram 
in one-loop. The determinant of A is given by 

detA = ^z 2 U (z x sx + z 2 s 2 + z 3 s 3 ){sl + s\ + s\ - 2sis 2 - 2s 2 s 3 - 2siS 3 ). 

(20) 

As an example, we try to calculate "Z exchange for tt vertex" i.e., 

2 2 2 2 

Pi = V 2 = mr, pi = s, 

mi = m 2 = WI4 = m5 = m, m 3 = rriQ = M, (21) 
m = 150(GeV), M = 91.17(GeV). 

Below the threshold, s < 4m 2 , the integral has no singularity and 
it can be. done easily by adaptive Monte-Carlo integration program 
BASEST Above threshold, the singularity, T> = , appears. It is 
hyperboloid of one sheet or that of two sheets in y space. The transition 
of the topology occurs at d = where d is defined by 

V = l gAg + c,' d = c- ^bA^b. (22) 



For the case in Eq. (|2T|) , the condition d — becomes 

z 3 = t, (M 2 -s)(s+3M 2 -Am 2 )t 2 +2M 2 (s-M 2 )t+M 2 (Am 2 -M 2 ) = 0. 

(23) 

By use of these formulas, we can calculate the value of integral based 
on the hybrid method as is described in detail in references. 




Figure 2: Cuts for non-planar vertex. 
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The calculation of the non-planar vertex at hand can be done easily 
by the dispersion integral. With respect to the variable s, we have a 
two-body cut as in FigJ^(a), and three-body cuts as in Fig.^(b) and the 
reversed one. 

j( s ) = I f ^ T(g/) , ds', ZT(s') = ST two + %T three (24) 
7T J s' — s — ie 

Here 

^sT two = J T (p 3 -> k\ + k 2 )dT 2 (k 1 , k 2 )T box (k 1 + k 2 V\ + P2) (25) 
and 

STthree = J T o(j>3 -> ki+k 3 +k 5 )dT 3 (k 1 , k 3 , k 5 )To(ki+k 3 +k 5 -> pi+p 2 ) + (reversed). 

(26) 

where tree amplitudes, T are expressed just by product of propaga- 
tor (s), and dT n stands for the n-body phase space. The one-loop box 
integral in Eq.(EH) is 



rte = J ' ixiviz \ 



-xyt — z(l — x — y — z)u + (x + y)M 2 + (1 — x — y) 2 m 2 ) 2 

(27) 

and is free from singularity because the denominator of is non-negative. 
Hence the 4-dimensional integrals in E q.(^5| ) and Eq. (|26|) are well be- 
haved. The singularity at s' — s in Eq. (|24j) can be handled by casting 
the integral into the form 

«J(.) = i-/ ar (' + g )- OT ('- ff W (28) 
2,71 J a 

By use of BASES, this is evaluated as 5-dimensional integral together 
all the integral variables, Feynman parameters and kinematical ones. 



Figure 3: Real part of non-planar vertex. 

In Fig.|3], we present the final results by these methods. The cal- 
culation of the integrals in Feynman parameter space is done at the 
same s repeatedly by changing the assignment of external momenta 
and masses in Eq. fl2"T|) . This diagnostics test works well for the points 
in the figures. 

In this report, we presented only one example for the real calcula- 
tion. However, the methods given here work in principle for any mass 
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parameters and external momenta" . The inclusion of numerator and 
the extension to four-point function are interesting and will be studied 
in the coming research. 
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llQne needs more sophisticated treatment if the infrared divergence exists. 



